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We propose to realize and observe Chern Kondo insulators in an optical superlattice with laser- 
assisted s and p orbital hybridization and synthetic gauge field, which can be engineered based on 
the recent cold atom experiments. Considering a double-well square optical lattice, the localized s 
orbitals are decoupled from itinerant p bands and are driven into a Mott insulator due to strong 
Hubbard interaction. Raman laser beams are then applied to induce tunnelings between s and p 
orbitals, and generate a staggered flux simultaneously. Due to the strong Hubbard interaction of 
s orbital states, we predict the existence of a critical Raman laser-assisted coupling, beyond which 
the Kondo screening is achieved and then a fully gapped Chern Kondo phase emerges, with the 
topology characterized by integer Chern numbers. Being a strongly correlated topological state, the 
Chern Kondo phase is different from the single-particle quantum anomalous Hall state, and can be 
identified by measuring the band topology and double occupancy of s orbitals. The experimental 
realization and detection of the predicted Chern Kondo insulator are also proposed. 

PACS numbers: 


Introduction - The search for new quantum states with 
nontrivial topology has become an exciting pursuit in 
condensed matter physics. The recent important ex¬ 
amples are the theoretical prediction and experimen¬ 
tal discovery of the 2D and 3D topological insulators 
(TIs) [I|,i. Recently, a few heavy-fermion materials (e.g. 
SmBe) were predicted to be candidates of time-reversal- 
invariant topological Kondo insulator (TKI), which orig¬ 
inates from the hybridization between itinerant conduc¬ 
tion bands and strong correlated / electrons HH. It 
was shown that screening the local moments by conduc¬ 
tion electrons leads to an insulating gap in the Kondo 
regime. The proposed scenario of TKIs is consistent with 
the transport measurement , angle-resolved photoe¬ 
mission spectroscopy M and scanning tunneling spec¬ 
troscopy [I10. Note that the essential difference be¬ 
tween a TKI and conventional TIs is not topological clas¬ 
sification, but the strong correlation physics for / elec¬ 
trons which are absent in conventional TIs. To measure 
such strong correlation physics can directly distinguish a 
TKI from conventional TIs, while this might be a chal¬ 
lenging task for condensed matter systems. 

In this letter, we propose to realize and observe 
a strongly correlated quantum anomalous Hall (QAH) 
phase, called Chern Kondo (CK) insulator, in an opti¬ 
cal lattice, motivated by the recent rapidly developing 
new technologies for cold atoms. These technologies in¬ 
clude the Raman coup ling scheme used to create 

spin-orbit interactions 17-E^. laser-assisted tunneling or 


shaken lattices to generate synthetic gauge fields 
241, and optical control of Feshbach resonance 2H-H|. 


*Materials and Correspondence should be addressed to: xiongjun- 
liu@pku.edu.cn 


Compared with solid state systems, the cold atoms can 
offer extremely clean platforms with full controllability 
to study many-body physics and topological phases |3^ 
1^ . Here, we consider a double-well square lattice with 
Raman-coupling-assisted s-p orbital hybridization to ob¬ 
serve CK insulating phases. Due to the strong Hubbard 
interaction of s orbital states, the Kondo screening is 
achieved when the applied Raman coupling exceeds a 
critical value, and then a nonzero renormalized s-p or¬ 
bital hybridization drives the system into a fully gapped 
CK phase. The parent state of s orbital is a Mott insula¬ 
tor which has no double occupancy due to the Hubbard 
gap, and this property keeps in the CK regime, signifying 
an essential difference from single-particle QAH insula¬ 
tors. With this unique feature the CK insulators can be 
identified by measuring band topology and double occu¬ 
pancy of s orbital in experiment. 

Model—F ot the realization we introduce a novel bi¬ 
partite square optical lattice created with the setup in 
Fig. [IKa), which has a staggered energy difference be¬ 
tween the neighboring A and B sites. This double-well 
checkerboard lattice can be realized based on NIST ex¬ 
periments [ 4 ^, |4l|, and more details will be presented 
later. The minimal model for realizing the CK insula¬ 
tors includes several basic ingredients. First, the sublat¬ 
tices A and B are anisotropic along the local coordinates 
X jY = ±x -b y. The relevant orbitals for our consid¬ 
eration are the s orbital at A sites and px orbital at B 
sites [Fig. [Hb)], with the intraorbital AA/BB hopping 
couplings along the diagonal directions. In the moderate 
to deep lattice regime we can verify that taso is much less 
than tppo-, leading to a relatively flat band for s-orbital. 
Accordingly, in such regime the transverse tunneling 
between jpy-orbitals is much weaker than axial tunnel¬ 
ing tp [ 43 , [ 4 ^ . Thus for the practical consideration, we 
neglect the 7r-bond hoppings of the s and px orbitals 
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FIG. 1: (color online) (a) Sketch of generating double-well 
checkerboard superlattice by an incident laser beam with in¬ 
plane (y) and out-of-plane (z) polarized components. The 
S-component is partially reflected by the mirror P. This con¬ 
figuration generates a 2D checkerboard lattice Vi^e,tt{x,y) = 
— Vbjcos^ffca:) -1- cos^ (ky)] — Vi sin^[|(a; — i/)] — V 2 cos^[|(a; -I- 
3/)] [5| - (b) The lattice potential profile for {Vb,Vi,V 2 } = 
{8, 8, 5}Eii, with Er being the recoil energy, (c) Two Raman 
lasers are applied to induce s-px orbital hybridization and 
create staggered flux pattern in (b). 

by setting = tf = 0, and for the tr-bond hopping 
we have \tp\ » |t^|. Secondly, similar as the recent 
experimental studies we consider an onsite en¬ 

ergy difference AE = Eb — Ea to suppresses the bare 
tunneling couplings between s and px orbitals. On the 
other hand, the neighboring s-px hybridization can be 
indnced by a standard two-photon Raman coupling in 
the following form Vr = Vm cos{Su!t + kay), where Vm is 
the coupling amplitude, Scu and k/j = kny represent the 
differences of frequencies and wave vectors between two 
Raman photons, respectively [Fig.[TKc)] pll-l2^. The Ra¬ 
man potential drives the neighboring hopping between s 
and px orbitals when Slo compensates their energy dif¬ 
ference AE. Moreover, this laser-assisted hopping along 
y direction is associated with a phase (po = k^a {—kna) 
if the hopping is toward (away from) B sites with A be¬ 
ing lattice constant, rendering a staggered flux pattern 
depicted in Fig. [T] (b) with the flux $ = 2(/)o. Finally, 
we consider a (pseudo)spin-1/2 two-copy version of the 
s-px band model, so the Hubbard interactions can be 
introduced for the spinful fermions (e.g. ®Li or ^°K) and 
denoted by Us (Up) for the s {px) orbitals. With these 
ingredients we can write down the effective Hamiltonian 
hy H = Hq + Hint, where 

^0 = ^\a^i±Ya ~ ^ss\„Sia + Pxio-TxiibAcr] 

IG 

+ XI + h.C.] (1) 

-^int — Us ^ ^ Up ^ ^ ^pxj't^Pxjj-' (^) 

* j 

Here = AE — Suj is the two-photon detuning for the 
Raman coupling [Fig.[TKc)], and pxja {sla ^^^Pxja) 


are annihilation (creation) operators of s and px orbitals, 
respectively, fisia = and fip^js^ = P^xjoPXoo The 

Raman coupling induced hybridization satisfies E{y) = 
±tspU‘'^° for r = Eaij and F(r) = Etgp for r = ±ax, with 
tsp being the induced hybridization strength. 

In the single particle regime the present setup can re¬ 
alize QAH effect [s^. We rewrite Hq in the k space 

Ho = with Ckcr = isua,PxkaV- The 

Bloch Hamiltonian takes the form 


'Ho(k) = do(k)To -b da:{k)Ta; + dy(k)Ty + 4(k)T^, (3) 


where dx(k) = 2tspSinkySin(j)o,dy{k) = 2tspSuikx + 
2tsp sin ky cos (po, do/^Ck) = coskx — tj cos ky — 
As/2, and the Pauli matrices Tx,y,z act on orbital space. 
The single-particle spectra, given by E^ = (io(k) ± 
tjdxiky + dyikY -b dziky, has node points if c/q = rmr 
with m £ h, and is gapped when </> yb rrnr and yb 
±2(t^ -b t^). The gapped single-particle system is in 
the QAH phase if |As| < 2(t^ -b tj) Q. On the other 
hand, as studied below, a positive detuning A^ > 0 will 
be applied to reach the CK phase. 

Strongly interacting regime and Kondo transition.-To 
study the Kondo phase, we consider the strongly in¬ 
teracting regime for s-orbital states with Us/t^ ^ 1, 
and a weak interaction for px orbital states, satisfying 
Us/tY Up/tp . To treat the interaction appropriately, 
the weak Up terms are decoupled by mean-field approxi¬ 
mation flp^j-tflp^ji ^ npxjaflpxja IT-pxjt^Pxji- 1^1 
the strong interacting limit Us/t^ 1, the state of 
double-occupation on s orbital is prohibited. The 
system can be studied by the slave boson mean field 
theory with vanishing doublon [i^. We introduce the 
slave-boson operator 6j(6i) at i-th site to describe the 
creation(annihilation) holon state °s. Moreover, we con¬ 
sider Zcrsl^pxja instead of sl^pxja in eq. (U) with = 


bB„ 


— , ^ _ Here the slave operators in uniform 

vTTir^l_h2_j| 

mean field approximation are replaced by time- and site- 
independent numbers, and b. The hopping tj is also 
renormalized by a factor of Z^, which preserves Lut- 
tinger volume and captures the correct band renormaliza¬ 
tion. To eliminate the redundant many-particle config¬ 
urations of s orbital, the constraints A + 6^ — l) 

J2ia- {fisia - sl) with A, A^r being the Lagrange 
multipliers are introduced. The functional for free en¬ 
ergy per site takes the form M | 


^ ko-,± 



where /3 ^ = k^T, the quasiparticle spectra of the 
lower (—) and upper (-b) hybridized bands are given by 
“b £sk(T T \/ i^px kCT ^sk^r)^ ~b 4| Z(yi2]^ p J /2, 
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respectively, with ep^ko- = 2t^ coskx + Upnp^g, 

Esko- = -2cos ky - As + Xa and Vk = 

2itsp (sin exp [z^q] + sin ). The minimiza¬ 


tion of the free energy with respect to the mean field pa¬ 
rameters are determined by ^ = 0 with ^ = Sa,b, Xa, X, 
which yields several coupled saddle-point equations 


J 


r, . -ST^rZadZa ZadZa^ 1 

{A-A,.,A} = y;{ }x-y: 


COS ky T 




= —y 

2N ^ 


k,± 


It 


k,± L 

‘^pxkfT ^sk(T 


•\/(£pxkcr ~ £skcr)^ T dl^crl^kP 

I- 


(^pxkcr ^skCT)^s COS ky T |Tk| 

\/(Spxkcr — Sska)^ T 4|^crVkp 


fiEt) 


( 5 ) 
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FIG. 2: (color online) (a) Kondo temperature Tk = , direct 

band gap A and Chern number Chi as a function of hopping 
integral tsp with {tp , As, Up) = (1,0.1, 2, 2) in unit of tp 
for gauge field = 0 and Band dispersion with 
tsp = 0.8tp , cj)o = 0 (b) and cfio = ^ (c). Each band has two 
fold degeneracy with respect of spin projection a. 


where f{E) = T 1] is the Fermi-Dirac dis¬ 

tribution function with /i the chemical potential. The 
half-filling condition for the entire system is fixed by the 
equation: 1 = [fiEka) + /(^^ka)] /‘^E, which sets 
the chemical potential in the self-consistent calculation. 

The Kondo transition can be studied in the regime 
with As > 0, which implies that the onsite energy of the 
s-orbital is lower than the Fermi energy of the half-filled 
pjc-band. The direct hopping from px orbital to s or¬ 
bital is however forbidden due to the Hubbard gap in the 
s states, while the process of “cotunneling” which is a 
second-order virtual process can occur between s and px 
orbitals. This process is known as the Kondo coupling. 
In the weak Raman coupling regime, the s-orbital states 
stay in a Mott insulator decoupled from the itinerant 
Px band, since the spectral weight of s-band vanishes at 
the Fermi energy. We thus expect that the emergence of 
the mean-field effective hybridization, namely, the Kondo 
coupling between s and px orbitals requires that the 
laser-assisted tunneling tsp exceeds a critical value, be¬ 


yond which the Kondo screening is achieved and a fully 
gapped CK phase emerges. 

We show the numerical results in Fig. [2] by solving 
the mean-field equations ([SI) self consistently. We can see 
that the effective hybridization between s and px orbitals 
is strongly renormalized by a factor Za due to the Hub¬ 
bard interaction in the s orbital states. A critical laser- 
assisted tunneling is observed and when tgp < t^p, 
the effective s-px hybridization renormalizes to zero [see 
the red curves in Fig. [2)[a)]. This is in sharp contrast 
to a single-particle system. Accordingly, increasing tsp 
to exceed t^p, the coherent hybridization is gradually de¬ 
veloped, leading to the formation of heavy quasi-particle 
bands as shown in Fig.[^(b) and (c). The local moments 
are screened by spin flipping through the cotunneling pro¬ 
cess. A direct bulk gap Agap > 0 opens between the 
heavy quasi-particle bands for ((jq = 7 r/2. In contrary, if 
the gauge flux (j> = 0, the da:(k) in eq. (O vanishes and 
the bulk has gapless nodal points. 

To analyze the topology of the bands, the Chern num¬ 
ber for the n-th band is explicitly calculated by 

Chi = f {fiEk.)-fiEi:^a)] 

Jbz ,, s 

m{^n),a 

1 Im 


where Vx/y = d'Hyldkx/y are velocity operators along 
X jY directions. Fig. shows nonzero Chern num¬ 
ber for (j)Q = 7r/2. The nontrivial topological property 
originates from the band inversion between the s and px 
bands. The phase diagrams versus gauge field (p and tsp 
are shown in Fig. [3] (a,b), with the subfigure showing the 
gauge field dependence of Tk and Chern number at a line 
cut of tsp = O.Stp . The Tk and the direct gap Agap show 
symmetric symmetric behavior with respect to pQ = 0, 
while Chern number Chi changes sign across this point. 

Experimental realization and detection.-Now we turn 
to the model realization and detection of the CK insulat¬ 
ing phase. The square superlattice potential lAatt(a:, 2 /) 


[{n\^a\vx |TOk(T)(mkcr|'()y [nka)] 


( 6 ) 
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FIG. 3: (color online) Main panel: (a) Image map of Kondo 
temperatnre Tk oc IP and (b) direct band gap A as a 
function of hopping integral tsp and gauge field (j>o with 
{tp , As, Up) — (1,0.1, 2, 2). Subfigue: (a) Kondo temper¬ 
ature Tif (open pentagon) and (b) band gap A (open square), 
Chern number Chi (open circle) versus (jio at tsp — O.Stp . 


can be generated through the setup shown in Fig. [II (a ) 
with an incident laser field E = (Eiy + £’ 2 z)e®^““^l[. 
Neglecting all the irrelevant phase factors, the in-plane 
polarized (x,y) components generate the standing wave 
as Eixy = 2 £’i [cos(fcoa^)y + cos(fcoj/)^] • The out-of¬ 
plane polarized (z) component generates the light field 
as E^z = zE 2 (E'"°^ + + ^^-ikox+z^ ^^xkoy+iny 

Here we set that the z component light is partially re¬ 
flected with a ratio a by the mirror P. The tt phase 
is the difference of phases acquired by x and z compo¬ 
nents of the laser beam traveling from the lattice center 
to mirrors M 3 and P, respectively, and back to the lat¬ 
tice. From Exy and E^z one can obtain the total lattice 
potential Vi,B±t(x,y) with the amplitudes {Vo,Vi,V 2 } oc 
{djEip — 4a|E2p,4(I -I- a^)|i52p,8a|E2p}. Hence the 
amplitudes Vop ,2 can be readily controlled by tuning a 
and Ei^ 2 - For example, taking that a = 0.35 and \Ei\'^ = 
1 . 47 |i? 2 p, we have Vq = SEr and Vi = 9,Er if setting 
V 2 = ^Er. This regime gives tp « 0.1ER,t^ k. O.OIE/}, 
and AE = Ep^ — Eg fa 1.75 Er, which well meet the 
former two ingredients of realization. 

The configuration for a strong interaction Ug at A sites 
and weak Up at B sites can be reached with optical con¬ 
trol of Feshbach resonance, which manipulates the atomic 
interactions by optical field strength 2^31|. The stag¬ 
gered modulation of interactions at A and B sites can be 
tuned with a periodic optical potential which minimizes 
(maximizes) at A (B) sites. Interestingly, in our case the 
Hi_ 2 -terms in the optical lattice potential match the re¬ 
quired profile and can play such role. We note that such 
a spatial modulation of interactions has been achieved in 
the recent experiment with a long life time 311. 

The full controllability of cold atoms can enable us to 
distinguish the CK insulating phase from single-particle 
QAH states. The present CK insulator can be identified 
by three characteristic features, namely, a critical s-px 
coupling strength tgp for CK phase transition, the non¬ 
trivial topology in the bulk band, and the Mott behavior 


of s-orbital. The former two features can be confirmed 
by band topology measurement through Hall effect 
which exists for tgp > t^p but is absent for tgp < f^p, while 
the last one is a key feature which makes the CK phase 
be exceptional. Since the parent state of s-orbital is a 
Mott insulator, the double occupancy of s-orbital is sup¬ 
pressed by the Hubbard gap. This property keeps both 
before and after the CK phase transition. Therefore, the 
third feature of CK phase can be detected by measuring 
the double occupancy of the s and px orbitals, which can 
be carried out in standard cold experiments 4^, 471. 


The double occupancy of s-orbital can be derived in 
the atomic limit. The partition function reads Z = 
2 " exp(—/3iJ„), where z = exp{/3p) and is the 
energy of n atoms occupying the s-orbital at an A site. 
We have that {Ei^,Ei,E 2 } = {0, —Ag,Ug — 2As} and 
the double occupancy of s-orbital Dg = z^ exp(—(3E2)/Z. 
Note that p, should be determined self-consistently with 
the total number of atoms in s and px states being fixed. 
On the other hand, the double occupancy of px-orbital 
can be obtained by Dp = (2/N)dP(Up,T,tgp)/dUp. 
From the numerical results shown in Fig [5] (a), one can 
find that Dg is greatly suppressed both before and after 
the CK phase transition, while Dp is unchanged versus 
tgp in the decoupled phase and increases in the CK phase. 
The later property reflects that the emergent hybridiza¬ 
tion in the CK phase leads to a net pumping of atoms 
from s-orbital to px orbital. For a comparison, in Fig [4] 
(b) we show the results for a single-particle QAH phase. 
With the same single-particle parameters, the double oc¬ 
cupancy Dg is much larger and exhibits qualitatively dif¬ 
ferent behavior upon increasing tgp. 

In conclusion we have proposed to realize and ob¬ 
serve CK insulators in an optical superlattice with laser- 
assisted s and p orbital hybridization and synthetic gauge 
field, which can be engineered based on the recently fast 
developing techniques in cold atoms. The predicted CK 
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FIG. 4: (color online), a) Double occupancy of px orbitals Dp 
and s orbitals Ds with (tp ,tj, As,Up,Us) = (1,0.1,2,0,20) 
as a function of tsp (in unit of tp) for various sets of temper¬ 
ature ksT. The circles at Dp curves represent the GK phase 
transition points, (b) Double occupancy Dp and Ds in single 
particle QAH phase with A^ = 0 and 2. 
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insulator can be identified by three characteristic fea¬ 
tures, namely, the existence of a critical s-p coupling 
strength for CK phase transition, the nontrivial topology 
in the bulk band, and the Mott behavior of the s-orbital. 
These features distinguish the strongly correlated CK in¬ 
sulating phase from the single-particle QAH states. 
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Supplementary Material — Chern Kondo Insulator in an Optical Lattice 


Double-well optical superlattice 


We provide here more details for the generation of the square superlattice potential Watt (a;,?/) through the setup 
shown in Fig. 1 (a) of the main text, with an incident laser field E = {Eiy-\-E 2 z)e^^°^. The initial phases are irrelevant 
and have been neglected. The light beam is reflected by the mirrors, with the forward and backward traveling fields 
intersecting and forming lattice in the middle region. In particular, the in-plane polarized (x, y) components are 
reflected by mirrors Mi, M 2 , and M 3 , and can form the standing wave in the form 

[ cos(fcoa: - tpi - V 2 )y + cos{koy + (f 2 )x] . (7) 

Here (pi represents the phase acquired by the in-plane light beam through the path from lattice center (intersecting 
point) to the mirror Mi, then to M 2 , and finally back to the lattice center, and (p 2 represents the phase acquired 
through the path from lattice center to the mirror M 3 [Fig. 1 (a)]. On the other hand, after traveling via Mi and 
M 2 , the out-of-plane polarized (z) component is reflected by a PBS (polarization beam splitter) and then partially 
reflected by P with a reflection ratio a. By controlling the optical path, we can set that the ^-component beam 
acquires an additional tt/2 phase relative to ip 2 by the in-plane polarized component traveling through the path from 
lattice center to PBS and then to P. This gives the standing wave for the z-component light 

_ 2 ^ 2 ( 6 *^““ + _|_ ^^-ikox+i2tpi+i2tp2+iTr _|_ ^^ikoy+iipi+i2tp2+i'^ 

It is easy to see that the phase factors ipi ^2 can be absorbed by shifting the position kox —>• kox + ^ 1+^2 and 
koU koy — ip 2 - Then the light fields read 

Ejjy = 2Ei[cos{kox)y + cos{koy)x], (9) 

= zE 2 (e*'=“="-be-*'=“^-bae-*'=‘>"^+”-bae*'=“^+"). ( 10 ) 

The above formulae imply that the generated standing wave fields are stable against phase fluctuations, which at 
most give rise to global shifts of the lattice and do not affect the model realization. From the above results we can 
get light field strength that (neglecting the constant terms) 

= 4|Eip[cos^(fcox)-bcos^(feo 2 /)], (11) 

= 4 |E 2 n(l -b a^)sin^[y(a: - y)]-b 2 acos^[Y(x -b 2 /)] - 4 a|E 2 |^cos^(fcoa;) - 4 a|E 2 |^cos^(fcoy) ( 12 ) 

Note that the lattice potential is proportional to the light field strength, with the proportional factor being the 
dipole matrix element. Therefore, the total lattice potential takes the following form 

VLa.tt{x,y) = - Vo[cos^{kox) +cos^{koy)] - W sin^[y (a; - y)] - W cos^[y (a; -b y)], (13) 

where the amplitudes {Vqj WiW} oc {4|Eip — 4 Q!|E 2 p, 4 (l -b a^)|E 2 p, 8 Q;|E 2 p}. Hence the amplitudes Vo, 1,2 can 
be readily controlled by tuning a and Ei^ 2 - This is the lattice potential considered in the main text. Taking that 
a = 0.35 and |Eip = 1 . 47 |E 2 p, we have Vq = SEr and Vi = 8 Er if setting V 2 = 7)Er. This parameter regime gives 
the lattice potential file shown in Fig. 1(b), with a large anisotropy along X/Y direction. The relevant parameters are 
obtained as Ri 0.1ER,tY ~ O.OIE^, and AE = Ep^ — Eg ~ 1.75E_r. It is noteworthy that the present realization 
is of topological stability, namely, the phase fluctuations in the incident laser beam cannot affect the configuration of 
the present double square optical lattice. 
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Slave boson theory 


The local correlation for the mixed s and px orbital lattice takes usual on-site Hubbard interaction 






UsY. 


nsj^nsj].. 


When the interaction Up for px orbitals is weak compared with band width of the dispersive px orbitals, the Hartree 


approximation hpxi^fip^ii ~ Y,o 




(i^Upxii is adequate to capture local correlation effects. The 


strong interaction Ug is treated more seriously by introducing auxiliary slave particles to describe the many-body 
configurations due to the localized nature of s orbitals. In the slave-boson theory the spin-cr fermion operator is 
written as 


^icr 





cLdi) R 


where the boson operators describe the holon bi, singly occupied states Cia-, doublon di, and fta- is a fermion operator. 
The explicit form of renormalization operators and takes the form 

Ricr — dj^di ^ 

_ l_ 

Rja — ^ 2(7 ’ 

recovering the Gutzwiller approximation as a saddle point solution of Hubbard model. 

To eliminate unphysical states, the Lagrange multipliers A and Acr are then introduced to ensure the local constraint 
of the completeness of the enlarged Hilbert space 


Qi — b\bi -|- + <i[di — 1 

(7 


and two equivalent ways to count the local fermion occupancy for each spin projection a between original fermion 
and slave boson representations 


Qia — d^di — 0 . 


The partition function within path integral formulism over coherent states of Fermi and Bose fields has the following 
form 


Z = Jv[b, &t] V [d, rft] V [c, ct] V [/, /t] V [A, A,] e- fo ^ 


where p ^ = k-^T and the corresponding Lagrangian C is given by 


^ ^ “ 1 “ ^ ^ f icT^'^ ^ ^ Up I ^ ^ '^pxif7^px‘i<7 I ^ ^ ^ Usdj^di 

i i<7 i \ <7 / i 

+ 4a^i±Yaflf^±Ya “ Px^aP Xr±Xa\ + 4a flaP X 3 ^ ^3 ,^+r + ^-C.] 


<b>c 


I T ^iaQi 


In the uniform mean-field slave boson approximation with infinite on-site Coulomb interactions C/g —>■ oo, the dou¬ 
blon vanishes to avoid the unphysical divergence and the boson operators are replaced by time- and site-independent 
numbers b, and d = 0. The resulting renormalization factor has a concise expansion 


bcr, 


Zi(7 — 


- cl\/l- 6^ - c; 


^^2 
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The functional for free energy per site is easily obtained from the partition function IF = —k^aTlnZ/N after integration 
of Fermi fields (/, /^) and takes the following form 

^ kcr,± Vo- / <y 


Here, are the lower (—) and upper (+) hybridized bands for the quasi particle with spin projection cr, respectively: 


— 1 
^ktr — 2 


^pxkcr “t“ £sl^(T i 'x/(^pxkcr S^skcr)^ F — 2^^ COS kx EUp'l^pxtr^ COS A:y A^-t- 

Acr and Vk = 2itsp (sin exp [i^o] + sin ). The extremization of the free energy with respect to the mean 

field parameters are determined by ^ = 0 with ^ = Co-, b, Ao-, A, which yields several coupled saddle-point equations 


{A-A.., A} = 


ZadZ„ 


Ccr' dCa' b db 


Zrr dZ„ 1 

vE 


cos ky T 


ct = 


1 

m 


E 

k.± 

= i-E 


It 


k,± L 

S^pxkfT ^sk(T 


\/ (^PxkfT S^sko-)^ H“ 4|EcrVk|'^ 


(^PxkfT ^skfT)^s COS J^Y |Vk| 
\/(^Pxkcr £^skcr)^ “1” 4|Eo-Vk|^ 

fiEtl 




-2 


where f{E) is the Fermi-Dirac distribution function. 

Note that in the main text, the substituted notion ^ and (^fia-yfja^ ■sL) > keeping the formalism 

complete in essence, is adopted without expanding the main text to introduce the auxiliary number and operator 
fxa- 


Effective model 


The standard periodic Anderson model is shown to connect the Kondo lattice via the well known Schrieffer-Wolff 
transformation . Below we shall derive the effective Kondo Hamiltonian describing the essential Kondo screening 
at =0 and Up = 0 limit. The original Hamiltonian H = Hp + Hg + Hgp takes the form 

Ep = E ^pP^XicrPxi+Xa- 

fo- 

Eg — 'y ^ E Ug y ^ 

i(7 i 

Esp — y ^ j^A'(r)s|g.pxjo-^yi+r T h.c. . 

<b)o’ 


Atomic projectors operating on f-th site s orbital are then introduced as follow 


Holon subspace : 7^i,H = (1 - ^sif) (1 “ 

Single occupied subspace : Vi^a- = ^sio- (1 — figia) 

Doublon subspace : Vi^o = ngi^figii 

Using identity + Fix> = 1, the hybridization between s orbitals and px orbitals can be written in 

projected subspace 

i7.p = E[^i(k)Ti«2(k)Th.c.] 

k 


with 


77i(k) = (k) 'Pi^as\„'Pi^}i-^px\^a exp[-zk • f], 

772 (k) = 'Pi,-Ds\„'Pi,s-^px\^G exp[-ik • i]. 


%G 
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A little algebra yields several commutation relations ['Hi{k),Hs] = As'Hi(k), [7^2(k), iJs] = (Ag — Us)'H2{k). For 
the localized s orbitals embedded in a conduction px band, the present model can be mapped into a spin exchange 
Hamiltonian through the explicit canonical transformation e~^ |5Cll | . The effective Hamiltonian up to quadratic order 
in (k) has a perturbation expansion via the Baker-Campbell-Hausforff formula and takes the form 

Ffeff = e^He-^ = H+[S,H] + ^^ [S, [S, H]] + O (V^ (k)) 

At the first order in Vk, the Schrieffer-Wolff transformation is constructed to eliminate the hybridization term with 


= E 

k 


-Ul (k) 



'H 2 (k) - h.c. 


By collecting the next leading order (quadratic order), we end up with the exchange Hamiltonian at large Ug limit 
-ffex = “ A -U ^ exp[ik' -i-ik- i]p^^^^,apxksl asi + k <—^ k'} 

kk' ^ ^ ^ ^ / i 

Fourier transforming back to real space, we arrive at a concise exchange Hamiltonian 

rF (r') 6j>^i+r'Pxj^PXj'4^^^ 

(b) <b') 

where J fa 2A {u‘-a ) • summarize, the minimal model Ffeff = Hp + iJex + ■ ■ ■ captures the essential physics. 
The exchange Hamiltonian 'Hex describes that the itinerant px bands hybridize with localized s orbitals, giving rise 
to the effective Kondo coupling. The screening of localized spin states of s orbitals we studied in the main text is 
adiabatically connecting the = 0 and weak Up limit. On the other hand, we remark that, in the main text, the 
Kondo screening is robust against the perturbative interaction Up which can be tuned by optical control of Feshbach 
resonance and then be treated by Hartree approximation. The essential cotunneling process is explicitly captured by 
the exchange Hamiltonian Hex- 


Double occupancy in atomic limit 


At the half-filling condition, to create a double occupancy in the lattice costs the energy U for the repulsive 
Hubbard interaction. Therefore, the strong repulsive interaction U can greatly suppress the double occupation of 
lattice sites and lead to Mott insulating phase. For the fermions in an optical lattices, the techniques to measure 
the fraction of double occupied sites have been well developed and demonstrated to detect the Mott insulating phase 
successfully 46, 4^. With a large repulsion Us, the hopping process for s orbital is greatly suppressed. In this case, 
we consider the atomic limit to approximately calculate the double occupancy by neglecting the kinetic energy. The 
partition function for i-th lattice site via Boltzmann factor has the following form 


Z, = ^ z” exp 


where z = exp [/3/i] is the fugacity, p is the chemical potential, n G {0,1, 2} labels three possible occupation numbers 
of fermions and Ei^n is the corresponding energy. We take the uniform approximation and drop the lattice site 
dependence. The local energy level for s orbitals is = {0, —A^,— 2A5 -|- Ug}- The chemical potential is 

located by the occupation number Ug in s orbitals 

2z exp ipAg) -\- 2z^ exp {2j3Ag — pUg) 

1 -I- 2zexp (/3As) -I- z^exp (2,5As - /3C/s) ’ 

Before the Chern Kondo phase transition, the chemical potential p is fixed at half-filling for s-band itself and the 
particle number ris = 1. On the other hand, after the Chern Kondo phase transition, p should be determined by the 
self-consistent solutions. The double occupancy reads 

z^ exp [-PUg] 

1 -I- 2z -I- z^ exp [—^Ug] 


D 
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where 2 ; = 


exp[-/3C/«]-l 

2exp[-^(7s] 


At high temperatures, the double occupancy saturates at Dg 


(V2+3ras-V2^^)'‘ 

16n| 


While, the asymptotic expansion at low temperature limit takes the form Dg = exp [—pUg] showing the exponen¬ 
tial behavior with respect to inverse temperature /3 and Hubbard interaction Ug. The double occupancy is negligible 
small at large Ug relative to temperature which is consistent with the vanishing doublon assumption in slave 

boson calculations. 
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